Fidelity Decay for Phase Space Displacements 
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In this letter we analyse the behavior of fidelity decay under a very specific kind of perturbation: 
phase space displacements. Under these perturbations, systems will decay following the Lyapunov 
regime only. Others universal regimes discussed in the literature are not presented in this case; 
instead, for small values of the perturbation we observe quantum freeze of the fidelity. We also show 
that it is possible to connect this result with the incoherent neutron scattering problem 
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Fidelity decay involves the overlap of two initially 
identical wavefunctions evolving under slightly different 
Hamiltonians 0. This is also known as Loshmidt echo. 
The subject has received much attention over the past 
few years as a model useful for discussion of decoher- 
ence and/or quantum classical correspondence |^,|^,^|^. 
For a generic perturbation of the Hamiltonian, it was 
shown that the decay follows three different universal 
regimes the Gaussian regime, the Fermi Golden Rule 
(FGR) regime and the Lyapunov regime, followed by one 
strong non-generic regime. The Lyapunov regime is espe- 
cially interesting because it is a way to measure the classi- 
cal Lyapunov exponent of a system from a pure quantum 
measurement 0; and it is a universal environment- free 
coherence decay . 

In this letter we analyze the behavior of Loshmidt 
Echo under a non-generic perturbation of the Hamilto- 
nian, namely phase space displacements. These Losh- 
midt types are suggested by experiments involving neu- 
tron scattering (where the phase space displacement is a 
momentum boost 0) and to some extent by electronic 
transitions in molecules and solids, where the displace- 
ment is along the position axis, with very little other 
change in the potential |1| . The non-generic character of 
this kind of perturbation might not be expected to follow 
the usual sequence of regimes, and indeed we show the 
phase space displacement presents only one of the univer- 
sal regimes, the Lyapunov one. In this letter we will fo- 
cus on incoherent neutron scattering. In general any kind 
of experimental measurement of momentum-momentum 
time correlation, position-position time correlation or a 
combination of both can be viewed as a Loshmidt Echo 
under phase space displacements. 

Under the formalism as presented by Lovesey Ig (which 
is based on Van Hove's Hnear response theory [ifll), the 
incoherent part of differential cross section for neutron 
scattering can be calculated from the correlation func- 
tion. 
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position operators of the nuclei and H is the typical 
Hamiltonian of the target system. In the Heisenberg rep- 
resentation, the correlation function can be written as: 



r,,(k,t) 



-jk-fj(0)^ik-fj(t) 



The ensemble average of the correlation function can be 
approximated by 0: 
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where \a) are coherent states in N degrees of freedom. 
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and <& (a) — e '^^fciC"), Introducing the 



notation e n^* = e *'<^'"je'?*e*'^''j, we can identified the 
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as the same kernel of a typical Loshmidt Echo problem: 
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where the brackets means ensemble average, fj are the 



Before we carry out a semiclassical evaluation of the 
fidelity decay in phase space displacement problems, it 
is interesting to analyze why it is possible to use the 
Diagonal Approximation as called in 0, and why this 
approximation is even better in the case analyzed in this 
letter. 

The Diagonal Approximation is based on the assump- 
tion that, given an perturbed Hamiltonian Hy = H + V, 
one can, classically, make the approximation that a tra- 
jectory 7 (i) in configuration space does not change, 
though its action is changed by AS"'' = Sy — S"' — 

dt'V [f {t')]. Unfortunately this assumption is not 
valid anymore for chaotic system due to the exponential 
divergence of single orbits. 

As pointed out by Cerruti and Tomsovic , and 
Vanicek and Heller |3] the reason this approximation 
is still valid for chaotic systems under a semiclassical 
formalism is mainly due to structural stability theo- 
rems Given a uniformly hyperbolic Hamiltonian 
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FIG. 1: Illustrative view of structural stability or orbits. Top 
plot: Structural stability for generic perturbations; 7 is a 
generic orbit of unperturbed Hamiltonian, 7y is the same ini- 
tial condition under perturbed Hamiltonian, 7' is a different 
orbit of the unperturbed Hamiltonian, and 'y'y is this different 
orbit under perturbation. Bottom plot: Phase space displace- 
ments; Labels are the same, except that subscript k is used 
instead of V for perturbed orbits. Notice that 7^ is barely 
visible since it is hidden under orbit 7. 



system, under a generic small perturbation it is almost al- 
ways possible to find a trajectory 7^ {t) of the perturbed 
system vifhich is arbitrarily close to the some unperturbed 
trajectory 7 [t) (even though they do not share common 
end points). Since semiclassical integrals involve a set 
of neighboring trajectories, both 7^ [t) and 7 (t) are in- 
cluded. In the classical limit of the Loshmidt Echo, this is 
equivalent to saying that the particle follows 7 (t) in the 
forward direction, but 7^ (t) in the backwards direction 
(i.e., that is the path which minimizes the action!). 
If the perturbation is given by a uniform boost in 



phase space, the diagonal approximation is exact apart 
from a time independent correction to the action: For 
any given trajectory 7 {t), with end points (qo, Qf ), there 
is a trajectory of the perturbed Hamiltonian H\^, 7^ it) 
with end points (qQ,qf) that follows exactly the previ- 
ous path, aside from small corrections at the extremities. 
The correction to the action is small and time indepen- 
dent. The fact that each individual pseudo-trajectory [l^ 
7 [t) -f 7y [t) has a time dependent correction to the ac- 
tion is fundamental to the existence of the FGR regime in 
generic Loshmidt Echo system. This is an intuitive rea- 
son why the FGR regime is absent in the case of phase 
space displacements. In Fig. ^ this discussion is illus- 
trated. 

For a more quantitative approach to the problem, we 
want to semiclassically evaluate 
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where \a) is a coherent Gaussian wavepacket, 
(r|a(ro,qo)) = (w2)~^e*P-("-'-») 



Following [13, we can approximate the semiclassical 
propagation of a coherent Gaussian wave packet using 
the Van Vleck propagator and expanding it around Vq: 
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where p = ^dS/dr\^^^^, the classical momentum for a 
trajectory starting on Tq and finishing on r' , Sj is the 
classical action along the same path, Vj the Maslov in- 
dex and Cj = \—d'^Sj (r, r'; t) /dridr'j |^_^ , all calculated 
over a given allowed classical trajectory 7j and summed 
over all possible classical trajectories. 

The kernel / (t) involves a double sum over all classical 
trajectories, j and j' , such that it can be interpreted as a 
overlap between two wavepackets, one that is boosted 
and subsequently propagated and another one that is 
propagated and subsequently boosted 0. However, as 
pointed out before, one can assume that only the overlap 
between the two trajectories that minimize the action of 
the pseudo-trajectory contributes significantly, while the 
others oscillate strongly. The two trajectories have the 
same action and the same Maslov index, so we can just 
assume that j = j' and write: 



e-'^-^^x 
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The last equation involves finding all possible classical 
paths which return to their initial position after time t. 
These can be large in number, and worseLthe prefac- 
tors Cj can diverge. Vanicek and Heller in |3| introduce 
the idea of using the IVR representation [lj| in a similar 
integral. Because there are actually two semiclassical in- 
tegrals simplified in the single integral above, instead of 
the usual term ^/Cj, we have Cj. However the term Cj 
is Jacobian of the transformation and disappear after the 
change from final position representation to initial mo- 
mentum representation. After transformation is done, 
we do not only eliminate the divergences in the integral, 
but we do not even need to calculate the term Cj any- 
more. We do not need also to consider more that one 
path, since the integral in initial momentum representa- 
tion has a unique allowed classical path. Introducing a 
new variable ^ = (p — Po) /ft— k/2, the integral becomes, 
after a change to initial momentum representation: 

I{t) = (^-^y e-'^ -e-^ J dV-'e"-^«^ 

where r' — r'(ro, p (po, ^) ;t) is the final position for a 
particle with initial position Tq and initial momentum p. 
M (t) is then given by: 




For a chaotic system, trajectories with a small initial sep- 
aration Aqo in phase space, i.e. Aqg ^ ALe^^ (where 
AL is the size of phase space and A the leading classical 
Lyapunov exponent of the system) we have that: 

|r; - r^l - Aq - ^l^^L^e^ = - 61 e^. 

On the other hand, if Aqo > ALe^^, then trajectories 
are uncorrelated and the mean separation is not depen- 
dent on time. In general, M (t) comes in two pieces, cor- 
related and uncorrelated. The uncorrelated part (t) 
can be calculated assuming a random phase approxima- 
tion in the integral above, 

Mu{t) = e 

that is, the uncorrelated part is a small time independent 
decay. In the generic Loshmidt Echo problem, the decay 
due to the uncorrelated trajectories is recognized as the 
responsible for the FGR part of the fidelity decay 
The result above means that there is no FGR regime in 
the fidelity decay of phase space displacements. 

Before we calculate the correlated part of the decay, we 
should notice that there is also no Gaussian decay in this 
problem. The Gaussian decay in generic Loshmidt Echo 



problems is obtained through first order quantum me- 
chanics perturbation theory. However, the displacements 
in phase space do not change the spectrum of the system 
(besides some possible shift of the ground state), so all 
the first order correction are zero, as the Gaussian contri- 
bution to the fidelity decay. This is actually analogous to 
the problem of spectrum change due to perturbation in 
the shape of billiards . When the billiards are simply 
displaced, there is decay of the overlap of eigenfunctions 
for example but no change in the spectrum. 

In order to calculate the correlated part of the fidelity 
decay, we will approximate the term exp («k • {r[ — r^) 
by its average over phase space. Because of the symmetry 
of phase space, r[ —v'2 will approach a Gaussian centered 
on zero with variance: 

((r;-4)^)~a4(4i-6)^^^*. 

So, 

and the correlated part of fidelity decay become: 
Mc{t)K{ — \ e ^ X 

«e-^(l + 2(ka)^e^^*)^ 
V2|k| cr' 

The last line reveals that the Lyapunov decay still exist- 
ing for the phase space displacement and it is actually 
the only time dependent fidelity decay for this kind of 
perturbation. The only condition we required upon the 
systems is the classical uniformly hyperbolicity. This is 
the first example of a pure Lyapunov decay of fidelity; 
it is also experimentally realizable. Because the fidelity 
decay M {t) is already averaged over coherent states, its 
thermal average is itself, {M{t)). So, it is straight for- 
ward to show that the autocorrelated function as defined 
in is just: 

r,,(k,t)«|(/(t))|-v/MlI). 

As a check and illustration of this work, we numeri- 
cally calculated the fidelity decay for a quantum stan- 
dard map with phase space displacement perturbations 
to compare with our predictions made with the semiclas- 
sical formalism. The results are shown in figure [3 for dif- 
ferent values of momentum boost k and different values of 
kick strength of the map Kg (which means different Lya- 
punov exponents). All numerical calculations were made 
with N = 5000 states. For small values of the "pertur- 
bation strength" (fc), we notice that the time-dependent 
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In between, there is a transition re- 



gion where the fidehty is frozen at intermediate values. 

In conclusion, we showed that interesting things hap- 
pen when we focus on particular form of perturbation 
in Loshmidt Echo problems. In particular, if we choose 
phase space displacement as particular form for the per- 
turbation, all the time dependent contribution to the fi- 
delity decay vanishes beside the Lyapunov contribution. 
This is not only an interesting example but also a nat- 
urally arising perturbation in the incoherent part of the 
differential cross section in neutron scattering. 

D. V. Bevilaqua was partially supported from CAPES, 
Coordenagao de Aperfeigoamento de Pessoal de Nivel Su- 
perior and NSF grant no. CHE-0073544. 




FIG. 2: Top Plot: Fidelity decay for different perturba- 
tion values. Dotted, dashed and dot-dashed lines are nu- 
merically calculated for the standard map with boost k = 



10-^10" 



10- 



respectively. Solid line is exponential decay 



at Lyapunov rate (where A — 2.7505) and "big"-dotted line 
is the "quantum" ergodic limit. All numerical calculations 
where done with kick strength Ks = 7.0, and averaged over 
20 coherent states. Bottom Plot: Fidelity decay for different 
kick strengths. Solid and dotted lines are Lyapunov decays 
with A = 2.7505 and 4.6497, respectively {Ks of 7.0 and 18.0). 
The dashed and dot-dashed lines are numerically calculated 
with Ks of 7.0 and 18.0, respectively, and both with k = 10~^. 
Inset on bottom plot shows the numerical calculated fidelity 
decay for k = 10~*, with dashed line correspoding to Ks — 7.0 
and dot-dashed corresponding to Ks = 18.0 



behavior of the decay is suppressed and the fidelity is 
frozen in a 'plateau' that is independent of the param- 
eters of the map (Ks) For bigger values of fc, the 
results seen to follow a Lyapunov decay until it reaches 
the quantum ergodic limit. To be more precise, given a 
quantization of the momentum and mean separation be- 
tween momentum levels Ap, if the perturbation k >> Ap 
then the Lyapunov decay is observed. On the other side, 
ii k « Ap then the fidelity is frozen on plateau with 
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